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Negative-parity baryon spectra in quenched anisotropic lattice QCD 
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We investigate tlie negative-parity baryon spectra in quenched lattice QCD. We employ the 
anisotropic lattice with standard Wilson gauge and 0{a) improved Wilson quark actions at three 
values of lattice spacings with renormalized anisotropy ^ = Oc/a-r — 4, where and a-r are spatial 
and temporal lattice spacings, respectively. The negative-parity baryons are measured with the 
parity projection. In particular, we pay much attention to the lowest SU{3) flavor-singlet negative- 
parity baryon, which is assigned as the A(1405) in the quark model. For the flavor octet and decuplet 
negative-parity baryons, the calculated masses are close to the experimental values of corresponding 
lowest-lying negative-parity baryons. In contrast, the flavor-singlet baryon is found to be about 1.7 
GeV, which is much heavier than the A(1405). Therefore, it is difficult to identify the A(1405) to be 
the flavor-singlet three-quark state, which seems to support an interesting picture of the penta-quark 
iudsqq) state or the NK molecule for the A(1405). 

PACS numbers: 12.38.Gc, 14.20.Gk, 14.20.Jn 



I. INTRODUCTION 

The lattice QCD simulation has become a powerful 
method to investigate hadron properties directly based 
on QCD. The spectroscopy of lowest-lying hadrons in the 
quenched approximation, i.e., without dynamical quark- 
loop effects, has been almost established, and reproduces 
their experimental values within 10 % deviations |lj. Ex- 
tensive simulations including dynamical quarks are in 
progress and would give us detailed understanding of the 
spectra of these ground-state hadrons 0]- In contrast, 
several lattice studies on the excited-state hadrons have 
been started 0, S IE IE S 0> very recently, and their 
calculations are far from established even at the quenched 
level. In this paper, using anisotropic lattice QCD, we 
investigate the low-lying negative-parity baryon spectra, 
particularly paying attention to the flavor-singlet baryon. 
Another purpose of this paper is to examine the corre- 
spondence between the flavor-singlet three-quark (3Q) 
state and the A(1405). 

In the context of the flavor-singlet baryon, the A(1405) 
is known to be one of the most mysterious hadrons. 
The A(1405) is the lightest negative-parity baryon, al- 
though it contains strangeness. Moreover, there are two 
physical interpretations on the A(1405). From the view- 
point of the quark model, the A(1405) is described as 
the flavor-singlet 3Q system. As another interpretation, 
the A(1405) is an interesting candidate of the hadronic 
molecule such as the NK bound state with a large bind- 
ing energy of about 30 MeV. We aim to clarify whether 
the A(1405) can be explained as a flavor-singlet baryon 
in quenched lattice QCD. 

Historically, excited-state baryons have been so far 



mainly investigated within the framework of the non- 
relativistic quark model, in which baryons can be classi- 
fied in terms of the spin- flavor SU{6) symmetry. While 
the ground-state baryons have completely symmetric 
spin-flavor wave functions and form the 56-dimensional 
representation, the low-lying negative-parity baryons are 
parts of the L = 1 orbitally excited states and belong 
to the SU{6) 70-dimensional representation. We sum- 
marize the classification of the SU (6) symmetry and its 
assignment to experimentally observed baryons in Table 
H Both the non-relativistic [l^ and the semi-relativistic 
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11] quark models reproduce the negative-parity baryon 
spectra fairly well in the octet and the decuplet sectors 
[l3 |. Such success of the quark model implies that the 
constituent quark picture holds well and the gluonic ex- 
citation modes play a less important role. The potential 
form is being clarified with the recent lattice QCD calcu- 
lations of the static 3Q potential fl^ IT^ . Furthermore, 
the large gluonic excitation energy |l5l | obtained with lat- 
tice QCD explains the reason why the quark potential 
model without such gluonic excitations well describes the 
hadron spectra. 

Among the low- lying negative-parity baryons, the 
A(1405) is an exception of such success of the quark po- 
tential model. In fact, the A(1405) is much lighter than 
the lowest-lying non-strange negative-parity baryons, 
iV(1520) with = 3/2- and N{1535) with = 1/2". 
There are two physical interpretations proposed for the 
A(1405): an SU{3) flavor-singlet 3Q state, and an NK 
bound state, i.e., a penta-quark (5Q) system. The simple 
quark model is based on the former picture, and it pre- 
dicts that the A(1405) and the A(1520) with = 3/2" 
are nearly degenerate (10^.11). In this picture, the large 
mass difference between them is explained to originate 
from a large LS'-force , but such a strong LS splitting 
is not observed in other baryon spectra, and therefore 
it seems difficult to reproduce the mass of the A(1405) 
within the simple quark model. 
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Another interesting interpretation for the A(1405) is 
the penta-quark (5Q ) sy stem or the NK bound state as 
a hadronic molecule HE Ell- Note here that the A(1405) 
lies about 30 MeV below the NK threshold, and this 
binding energy of about 30 MeV is rather large in com- 
parison with about 2.2 MeV, that of the deuteron. If 
this picture holds true, the large binding energy between 
N and K results in a significant role of the attractive 
effect for the K put inside nuclei or nuclear matter. In 
this way, the study of such an exotic and strange baryon, 
the A(1405), is important also for understanding of the 
manifestation of strangeness in the hyper-nuclei and the 
neutron stars. 

The 3Q and the 5Q states, however, would mix in the 
real world. Therefore, a more realistic question would 
be as follows. Which is the dominant component of the 
A(1405), the 3Q state or the 5Q state? We try to answer 
this question using lattice QCD simulations. In lattice 
QCD simulations, even if one chooses the operator as the 
3Q or the 5Q state, it generally overlaps with both the 
states through the quark and anti-quark pair creation. In 
the quenched simulation, however, owing to the absence 
of dynamical quark-loop effects, such a mixing between 
3Q and 5Q states are rather suppressed, which would 
enable us to investigate the properties of genuine 3Q and 
5Q states in a separate manner. 

In this paper, we focus on the 3Q state and investigate 
whether this picture can explain the mass of the A(1405). 
Apparent discrepancy with the experimentally observed 
mass implies that the penta-quark state gives significant 
contribution to the physical A(1405) state. In practice, 
lattice QCD results suffer from various systematic errors. 
It is therefore essential to compare the mass of the flavor- 
singlet 3Q state with other negative-parity baryon masses 
as well as with lowest-lying baryon masses. 

It is also important to understand the gross struc- 
ture of low-lying negative-parity baryon spectrum in re- 
lation to spontaneous chiral-symmetry breaking in QCD. 
If the chiral symmetry is restored, such as at high tem- 
perature and/or density, the masses of a baryon and 
its parity partner should be degenerate. Spontaneous 
chiral-symmetry breaking causes mass splitting between 
positive- and negative- parity baryons. It is important 
to study nonperturbatively negative-parity baryons in 
terms of the parity partners of the positive-parity ones. 

Since the negative-parity baryons have relatively large 
masses, their correlators rapidly decrease in the Eu- 
clidean temporal direction. As a technical improvement, 
we adopt an anisotropic lattice where the teniporal lat- 
tice spacing is finer than the spatial one, aa jig . With 
the high resolution in the temporal direction, we can fol- 
low the change of the correlator in detail and specify 
the relevant region for extraction of the mass. Thus ef- 
ficient measurements would be possible. This approach 
is efficient also for other correlators of heavy particles, 
such as the glueballs j^^. (The anisotropic lattice is ex- 
tremely powerful for the study of the finite temperature 
QCD [l3.l20ll2l| . where the temporal distance is severely 



limited in the imaginary-time formalism.) 

In this study, we adopt the standard Wilson plaquette 
gauge action and 0{a) improved Wilson quark action, 
for which the sizes of errors are rather well evaluated 
(221 i23] . The simulations are performed on the quenched 
anisotropic lattices with renormalized anisotropy ^ — 
aaja-r = 4 at three lattice spacings in the range of 
a^^ ~l-2 GeV. For these lattices, the quark parameters 
were tuned and the light hadron spectrum was calcu- 
lated in order to estimate the effects of uncertainties due 
to anisotropy on the spectrum p3 |. 

This paper is organized as follows. In Sectional we 
summarize the anisotropic lattice actions used in this 
study. We show the numerical results on the negative- 
parity baryon in Section IIIII and discuss their physical 
consequences in Section lTVI The last section is dedicated 
to the conclusion and perspective for further studies. 



II. ANISOTROPIC LATTICE 

We employ the standard Wilson plaquette gauge 
action and the 0{a) improved Wilson quark action 
on anisotropic lattices. We briefly summarize the 
anisotropic lattice action. The gauge field action takes 
the from 




I-iReTrC/,j(x) 



1 ~ -ReTr[/,4(a;) 



(1) 



with /3 = 2Nc/g'^. Here, U^^, denotes the parallel trans- 
port around a plaquette in ^-v plane, 

Uf,,{x) = U^(x)U,{x + fL)Ul{x + vMix). (2) 

The gluon field is represented with the link-variable as 
t/^ ~ exp(— zga^yl^). The bare anisotropy coincides 
with the renormalized anisotropy =^ aaja-r at the tree 
level. 

Note here that the bare anisotropy is no longer the 
same as ^ due to the quantum effect, and one needs to 
measure ^ through some physical observables for each in- 
put value of 7(5. Although ^ is in general a function of 
gauge and quark parameters, (/3,7g) and (k^^p) respec- 
tively, on a quenched lattice, the calibrations of the gauge 
and quark actions can be performed separately. For the 
gauge action of the form Klassen nonperturbatively 
obtained an expression of 70 in terms of /3 and ^ with the 
accuracy better than 1% using the Wilson loops |22|. We 
adopt the same lattice actions as those in Ref.'235 which 
made use of the Klassen's result for the gauge action and 
also performed sufficient analysis for the quark actions as 
mentioned below. 

For the Wilson type quark action, 0(a) improvement 
is significant in quantitative computation of hadron spec- 
trum. Among several types of the anisotropic lattice 
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TABLE I: Quark model assignments for experimentally observed baryons in terms of the spin-flavor SU(6) basis [l^ . 



SU{6) rep. 


5(7(3)/ rep. 




5 = 


5 = -1,/ = 


5 = -1, J = 1 


S = -2 


5 = -3 


56 (L = 0) 


^8 


1 + 

2 


iV(939) 


A(1116) 


E(1193) 


H(1318) 








3 + 
2 


A(1232) 




S(1385) 


H(1530) 


0(1672) 


70 (L = 1) 


^8 


1 - 

2 


iV(1535) 


A(1670) 


E(1620) 


H(?) 








3- 
2 


iV(1520) 


A(1690) 


E(1670) 


H(1820) 






*8 


1 - 
2 


iV(1650) 


A(1800) 


S(1750) 


H(?) 








3- 
2 


iV(1700) 


A(?) 


E(?) 


H(?) 








5- 
2 


iV(1675) 


A(1830) 


E(1775) 


H(?) 








1 - 
2 


A(1620) 




E(?) 


H(?) 


m 






3- 
2 


A(1700) 




m 


H(?) 


m 






1 - 
2 




A(1405) 












3- 
2 




A(1520) 









quark action, we use the form proposed in Refs. [20l 
l24j . As a merit of this form, the cahbration with a good 
precision is rather easy in the light quark mass region, 
since the quark mass dependence is expected to be small 
there, as was numerically shown in Ref. p^ . 
The quark action is written as 

5f = ^V^(a;)/^(a;,y)^(y), (3) 
x,y 

K{x,y) = 5.^^y 

{(l-74)C/4(x)<5,+4_^ + {l+j4)Ul{x - i)S,_J,^y} 

i 

-KcrCE ^ cri4Fi4Sx,y - rK^CB ^ crijFijSx,y, (4) 

i i>j 

where tp denotes the anticommuting quark field, and 
Kt the spatial and temporal hopping parameters, respec- 
tively, r the spatial Wilson parameter and ce, cb the 
clover coefficients. The field strength i^^^ is defined with 
the standard clover-leaf-type construction. In principle, 
for a given k^, the four parameters r, ce and 

Cb should be tuned so that Lorentz symmetry is sat- 
isfied up to discretization errors of O(a^). Following 
Refs. [23, 113, 13 1 the spatial Wilson parameter 

as r = 1/^ and the clover coefficients as the tadpole- 
improved tree-level values, namely, 

r = l/^, CE^l/ucyul, CB^l/ul- (5) 

To reduce large contribution from the tadpole diagram, 
the tadpole improvement ;25j is applied by rescaling 
the link variables as Ui{x) —^ Ui{x)/ua and 6^4 (x) — > 
U4^{x)/ut, with the mean- field values of the spatial and 
temporal link variables, Ua and Mr, respectively. This is 
equivalent to redefining the hopping parameters with the 
tadpole-improved ones (with tilde) through = ka/uc 
and jur- We define the anisotropy parameter ^e 

as 7f = K,rlK„. This parameter in the action is to be 



tuned nonperturbatively in the numerical simulation. It 
is convenient to define k as 

^ = ^-2(7F + 3r-4) = 2(mo7F + 4), (6) 

where mg is the bare quark mass in temporal lattice 
units. This n plays the same role as in the case of 
isotropic lattice, and is convenient to parameterize the 
quark mass together with the bare anisotropy ^e- 

The above action is constructed following the Fermilab 
approach (26| , which proposes to tune the bare anisotropy 
parameter so that the rest mass and the kinetic mass 
equal each other. In practice, hadronic states are conve- 
nient to carry out this program. In Ref. ^23], the bare 
anisotropy was tuned nonperturbatively using the rela- 
tivistic dispersion relation of the pseudoscalar and vector 
mesons. The main result of Ref. 123] is as follows: They 
tuned 7i? in the quark mass range from the strange to 
charm quark masses with the accuracy better than 1%, 
and found that the tuned bare anisotropy, , is well fit- 
ted to the linear form in m^, where = ~'*7^)/2'? 
is naively defined quark mass. Then 7^, in the massless 
limit was obtained within 2% error, i.e., 1% as a sta- 
tistical error and 1% as a systematic error in the form 
of fit in terms of rriq. Then, they computed the light 
hadron spectrum using the value of 7^ at the chiral limit, 
and observed the effect of uncertainty in 7^, on the spec- 
trum for physical quark masses of the 1%-level. In the 
present study, we also treat the same quark mass region 
and therefore adopt the value of 7^ in the chiral limit. 
The precision of 2% error in 7^, is sufficient for the present 
purpose. 

As was pointed out in Refs. H^ll^l, with the choice 
r = 1/^, the action Q leads to a smaller spatial Wilson 
term for a larger anisotropy ^. Since the negative-parity 
baryons measured in this paper are the lowest state of 
the parity projected baryon correlators, the statements 
in Ref. |23| for light hadrons also hold in our calculation. 
Even for the coarsest lattice in our calculation, the cutoff 
~ 4.0 GeV seems sufficiently large to avoid the arti- 
ficial excitation due to the doublers in the ground-state 
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signals. At least, the finest lattice with ^ ~ 8 GeV 
would be sufficiently large to avoid the doubler effect. 



III. NUMERICAL SIMULATION 
A. Lattice setup 

The numerical simulation is performed on the same 
lattices as in Ref. Here, we briefly summarize the 

fundamental parameters and physical quantities. We 
use the three anisotropic lattices with the renormalized 
anisotropy ^ = 4, at the quenched level. The statisti- 
cal uncertainties are, otherwise noted, estimated by the 
Jackknife method with appropriate binning. 

The lattice sizes and parameters in generating the 
gauge field configurations are listed in Table ^1 The 
spatial lattice scales roughly cover 1-2 GeV. The 
values of bare gluonic anisotropy are chosen accord- 
ing to the result by Klassen |22|. The uncertainty of 
his expression is of the order of 1%, and in the follow- 
ing analysis we do not include this uncertainty in the 
quoted statistical errors. The gauge configurations are 
separated by 2000 (1000) pseudo-heat-bath sweeps, af- 
ter 20000 (10000) thermahzation sweeps at P=5.95 and 
6.10 (5.75). The configurations are fixed to the Coulomb 
gauge, which is convenient in applying the smearing of 
hadron operators. 

The mean-field values of link variables are determined 
on the smaller lattices with half size in temporal extent 
and otherwise with the same parameters for /3 = 5.75 
and 5.95, while at f3 — 6.10 the lattice size is 16^ x 64. 
The mean-field values, Ua- and u^-, are obtained as the av- 
erages of the link variables in the Landau gauge, where 
the mean-field values are self-consistently used in the fix- 
ing condition j2nl |. In a study of hadron spectrum, it is 
convenient to define the lattice scales through a hadronic 
quantity. We determine through the K* meson mass, 
niK* = 893.9 MeV(isospin averaged). The procedure is 
the same as in Ref. ;,23J, while with larger statistics. The 
result is quoted in Table Ull as a~^{mK')- 

The quark parameters are listed in Table IIIII These 
hopping parameters roughly cover the quark masses 
niq ~ ms-2TOs. The numbers of configurations are larger 
than those of the hadronic spectroscopy in Ref. As 
already noted in the previous section, the values of 
and taken from the result of Ref. (23]. Although 

the uncertainty of 2% level is associated with the values 
of "fF, the quoted errors of hadron masses in the follow- 
ing analysis do not include this uncertainty. According 
to Ref. [23, this uncertainty in the physical masses of 
vector mesons and positive-parity baryons are at most of 
order of 1%. For the negative-parity baryon masses, the 
effect is expected to be similar amount and not signifi- 
cant compared with the present level of statistical error. 
Therefore, this effect can be negligible in our calculation. 

As described later, we first extrapolate the vector me- 
son mass linearly in the pseudoscalar meson mass squared 



to the point at which the ratio of these meson masses are 
equal to the physical value, mK*/mK- At this point, 
aforementioned lattice scale is determined. The physical 
(m, d) and s quark masses are determined through the tt 



and K meson masses, 



139.6 MeV and mx = 495.7 



MeV (isospin averaged), respectively. 



B. Baryon correlators 

We measure the correlators in pseudoscalar and vec- 
tor meson channels and octet (S and A types), decu- 
plet and singlet channels of SU(3) fiavor representation 
of baryons. As listed in Table IIVI we use the standard 
meson and baryon operators which have the same quan- 
tum numbers as the corresponding baryons and survive 
in the non-relativistic limit. It is known that there are 
mainly two ways to choose the baryon operator. One 
is the operator taken here, {q'^Cj5q)q, and the other is 
of the form {q^ Cq)'y5q. There are two reasons why we 
take the former. It is well-known that the former opera- 
tors strongly couple to the the ground-state baryons and 
reproduce experimental values well. Therefore, it is suit- 
able for investigation of the parity partner of the ground 
state baryons. Furthermore, the recent lattice calculation 
shows that these two operators give the similar results for 
the negative-parity baryon spectrum while the latter is 



'e-pa 
0. 



more noisy 

For baryons, two of three quark masses are taken to be 
the same value as specified by the hopping parameter ki, 
and the other quark mass is specified by K2- This cor- 
responds to taking the same value for u, d current quark 
masses as uiu = rrid = m„. Then, the baryon masses are 
expressed as the function of two masses mi and 7712, or 
equivalently of ki and K2, like mB{K,i, ^2). In the source 
operator, each quark field is smeared with the Gaussian 
function of width ~0.4 fm. 

At large t (and large Nt — t), the baryon correlators 
are represented as 



X 

= (1 + 74) [cs+ • + bcB- ■ e- 



(Nt-t)mg- 



+ (1-74) • e-(^'-*)'"«+ + 



cb- ■ e 



, (7) 



where h = +1 and —1 for the periodic and antiperiodic 
temporal boundary conditions for the quark fields. Since 
we adopt the standard Dirac representation for 7 matri- 
ces, the upper and lower two components correspond to 
the first and second contributions of Eq. |7J . 

Combining the parity-projected correlators under two 
boundary conditions, one can single out the positive- and 
negative-parity baryon states with corresponding masses 
mB+ and m^- , respectively, without contributions from 
the backward propagating parity partners. In practical 
simulation, however, we take a sufficient temporal extent 
so that we can observe an enough range of plateau in 
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TABLE II: Lattice parameters in the gluon sector. The scale a'^^{mK») is determined from the K* meson mass. The mean- 
field values are defined in the Landau gauge. The statistical uncertainty of Ut is less than the last digit. The details for these 
parameters are described in j23|. 

(3 jG size Ua- Mr a'^^{mK*) [GeV] 

5.75 3.072 12^ x 96 0.7620(2) 0.9871 1.034( 6) 

5.95 3.1586 16^ x 128 0.7917(1) 0.9891 1.499( 9) 

6.10 3.2108 20^ X 160 0.8059(1) 0.9901 1.871(14) 



TABLE III: Lattice parameters in the quark sector. The values of -yp and Kc are taken from Ref. [2^ . 



p 






^corif 


values of k 


5.75 


3.909 


0.12640(5) 


400 


0.1240, 0.1230, 0.1220, 0.1210 


5.95 


4.016 


0.12592(6) 


400 


0.1245, 0.1240, 0.1235, 0.1230 


6.10 


4.034 


0.12558(4) 


400 


0.1245, 0.1240, 0.1235, 0.1230 



effective mass plot for extraction of mass in each par- 
ity channel, and hence there is no advantage in comput- 
ing correlators under two boundary conditions except for 
the reduction of statistical fluctuation. We obtain the 
baryon correlators at /3 = 5.75 under two boundary con- 
ditions, and compare the statistical fluctuations in the 
parity-projected correlator and in not projected one. We 
conclude that it is not worth doubling the computational 
cost and hence adopt only the periodic boundary condi- 
tion hereafter. Instead, at each /3 we obtain the corre- 
lators with the source at t = Nt/2 in addition to ones 
with the source at t = 0, and average them. This is ef- 
ficient to reduce the statistical errors for limited number 
of configurations. 

C. lattice QCD results for hadron masses 

Fig. n] shows the effective mass plots for the baryon 
correlators at = 6.10. The effective mass is defined 
without considering the contribution of the associated 
parity partner propagating backward from the source at 

We observe that in the region where the effective mass 
exhibits a plateau, the contribution of parity partner is 
sufficiently small. In particular for the negative-parity 
baryon channels, fine temporal lattice spacing seems to 
be helpful to specify the region in which the ground state 
dominates. 

The meson correlator is fitted to the single h ype rbolic 
cosine form and analyzed independently of Ref. |23l | . The 
results are shown in Fig. [3 and consistent with Ref. [23| . 
For the baryons, we fit the data to a single exponential 
form. The result is listed in Tables IVl IVll and rTlTl 

Following Ref. , we extrapolate the hadron masses 
to the chiral limit in terms of the pseudoscalar meson 



mass squared, instead of 1/k. The assumed relation be- 
tween PS meson mass and quark mass is 

mp5(TOi, 7712) = S • (mi -I- m2), (9) 

then for the degenerate quark masses, mi = m2, rn?ps — 
2Bmi holds. Instead of m^ (i=l,2), one can extrapo- 
late other hadron masses in term olmps{mi, rrii)'^ to the 
chiral limit. 

In our calculation for baryons, two of quark masses 
are taken to be the same value, mi, and the other quark 
mass m2 is taken to be an independent value. Then, 
the baryon masses are expressed as the function of mi 
and m2 like mB(mi,m2), and therefore they are to be 
depicted on the {1711,1712) plane. However, the result for 
the baryon masses seem to be well described with the 
linear relation, 

ms(mi,m2,m3) = mB(0, 0, 0) + • (2mi +m2). (10) 

Therefore, we fit the baryon mass data to the linear form 
in the sum of corresponding PS meson masses squared. 
The vector meson is also fitted to a linear function in 
mi -I- m2. 

Note here that, in quenched QCD, a non-analyticity 
appears in the chiral extrapolation near the chiral 
limit ■ For nucleons, it is reported that the depar- 
ture from the simple chiral extrapolation is observed for 
nijr < 400 MeV |23|. Although we have to keep this ef- 
fect in our mind, we do not argue it here because there is 
no distinct behavior for both the positive- and negative- 
parity baryons in our calculation with m^r > 600MeV. 

The results of fits for baryons are shown in Fig. |3| for 
each lattice. The horizontal axis is the averaged pseu- 
doscalar meson mass squared. 




(11) 



6 




o 



degenerate quark mass 



a 



too, 



16 




0.32 



0.30 



SS 0.28 Ozszszsa. 



S 0.26 



0.24 



0.22 



p=6.10 singlet(-) 
degenerate quark mass 



vnn. 



Onni 



360 




16 



24 



32 



40 



48 



56 



FIG. 1: Effective mass plots for octet, decuplet and singlet baryon correlators with degenerate quark masses at /? = 6.10. The 
symbols correspond to k =0.1230, 0.1235, 0.1240 and 0.1245 from top to bottom in each figure. The left figures are for the 
positive-parity baryons and the right for the negative-parity ones. 
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TABLE IV: Typical interpolating operators for various hadrons. For baryons, the contraction with the color index is omitted. 
C denotes the charge conjugate matrix. 



Meson 


Pseudoscalar 






M{K) = S75U 






Vector 






Mk{K*)^s-fkU 




Baryon 


Octet 




b^{tP) = {c^5)m 


[ua{disa^ — Sfjd-f) — 


daispU-, — U13S-,)] 




Octet (A) 


B.(A) = 




spd^i) + daispUj ~ 


U/SSj) — 2Sa{ui3d^ - dpUj)] 




Singlet 


Bc(Ai) 


= {C'y5)p-t[Uaidi3S^ - 


- spd^) + daispu^ - 


- U/sSj) + Saiupdj — d/sUj)] 




Docuplet 




1 = {Cjk)l3j[Ua{dl3Sj 


+ 8/3^7) + da{Sf3U~, 


+ U/3Sj) + Saiupdj + di3Ury)] 



TABLE V; Baryon spectrum at /3 = 5.75 in the temporal lattice unit. When quark masses are degenerate as ki = K2, the 
S-type and the A-type octet baryon correlators become identical. 









Positive-parity baryons 






Negative-parity baryons 








"loct{E) 


'Ttoct(A) 


TTising 






'T^oct{A) 




TTldcc 


0.1210 


0.1210 


0.4281( 9) 


0.4281( 9) 


0.6426(56) 


0.4606(16) 


0.5443(52) 


0.5443(52) 


0.5355(30) 


0.5630(54) 


0.1210 


0.1220 


0.4171(10) 


0.4186(10) 


0.6349(59) 


0.4518(17) 


0.5355(56) 


0.5356(57) 


0.5265(32) 


0.5548(59) 


0.1210 


0.1230 


0.4059(10) 


0.4091(10) 


0.6267(64) 


0.4433(19) 


0.5271(61) 


0.5270(64) 


0.5175(34) 


0.5467(66) 


0.1210 


0.1240 


0.3946(11) 


0.3998(11) 


0.6184(69) 


0.4352(21) 


0.5189(68) 


0.5178(78) 


0.5089(37) 


0.5386(78) 


0.1220 


0.1210 


0.4083(10) 


0.4066(10) 


0.6275(63) 


0.4430(19) 


0.5269(61) 


0.5269(60) 


0.5173(33) 


0.5465(65) 


0.1220 


0.1220 


0.3970(11) 


0.3970(11) 


0.6206(68) 


0.4341(20) 


0.5181(66) 


0.5181(66) 


0.5080(35) 


0.5382(71) 


0.1220 


0.1230 


0.3855(11) 


0.3874(11) 


0.6132(74) 


0.4256(22) 


0.5095(73) 


0.5091(75) 


0.4989(38) 


0.5300(79) 


0.1220 


0.1240 


0.3737(12) 


0.3780(12) 


0.6055(81) 


0.4176(25) 


0.5012(82) 


0.4994(90) 


0.4899(41) 


0.5217(93) 


0.1230 


0.1210 


0.3885(11) 


0.3845(11) 


0.6116(74) 


0.4259(22) 


0.5091(77) 


0.5099(74) 


0.4990(38) 


0.5302(81) 


0.1230 


0.1220 


0.3769(12) 


0.3747(12) 


0.6063(81) 


0.4171(24) 


0.5001(83) 


0.5008(81) 


0.4895(41) 


0.5218(89) 


0.1230 


0.1230 


0.3650(12) 


0.3650(12) 


0.6004(89) 


0.4085(27) 


0.4913(92) 


0.4913(92) 


0.4801(44) 


0.513(10) 


0.1230 


0.1240 


0.3527(13) 


0.3554(13) 


0.5939(99) 


0.4004(32) 


0.483(11) 


0.481(11) 


0.4707(48) 


0.505(12) 


0.1240 


0.1210 


0.3689(13) 


0.3612(13) 


0.5936(89) 


0.4099(30) 


0.489(11) 


0.4928(98) 


0.4805(46) 


0.514(11) 


0.1240 


0.1220 


0.3569(14) 


0.3512(14) 


0.591(10) 


0.4010(33) 


0.480(12) 


0.483(11) 


0.4708(49) 


0.505(12) 


0.1240 


0.1230 


0.3445(15) 


0.3414(15) 


0.588(11) 


0.3922(37) 


0.470(13) 


0.473(12) 


0.4610(53) 


0.496(14) 


0.1240 


0.1240 


0.3317(16) 


0.3317(16) 


0.585(13) 


0.3838(45) 


0.460(15) 


0.460(15) 


0.4509(58) 


0.486(16) 


fit 1 


range 


24-40 


24-40 


12-20 


28-40 


20-32 


20-32 


16-24 


20-32 



with Nq — 3 for baryons. The results of fits are displayed 
as the solid lines in these figures. The linear relation 
seems to hold well. 

As stated in Section ITTll we determine the scale 
through the K* meson mass. The physical (m, d) and 
s quark masses are determined with the tt and K me- 
son masses. The scales of the vertical and horizontal 
axes in Fig. |2| are set in this way. The corresponding 
hadron masses for the physical quark masses are listed 
in Table IVlTll These results for the meson masses and 
positive-parity baryon masses are consistent with those 
obtained in Ref. l23l. 



D. Systematic errors 

Finally before discussing physical implications of our 
numerical results, we briefly comment on the systematic 
uncertainties. 



• Anisotropy (calibration) 

According to the detailed inspection given in Ref. 
p^ . the 2% uncertainty in 7f causes uncertainties 
in hadron masses at 1% level. Although the effect 
of uncertainty coming from anisotropy on the neg- 
ative parity baryon masses is unknown, its size is 
expected to be the same level as for the positive- 
parity baryons and smaller than their statistical 
errors. Therefore, we do not perform a detailed 
analysis of the calibration uncertainty here. The 
uncertainty in 7^, the gluonic anisotropy parame- 
ter, is also kept within 1% level, and hence for the 
same reason as for jp, we do not argue its effect on 
the negative-parity baryon masses. 

• Finite volume effects 

Since the excited baryons may have larger spa- 
tial extent than the ground state baryons, they 
may seriously suffer from the finite size effects. 
Our present three lattices, however, have almost 
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TABLE VI: The same 


results as Table [V] for /3 — 


5.95. 












Positive-parity baryons 

'^oct(A) nosing 


notice 




Negative-parity baryons 

'T^oct{A) Tlsing 


TTldcc 


0.1230 
0.1230 
0.1230 
0.1230 


0.1230 
0.1235 
0.1240 
0.1245 


0.2785( 7) 
0.2722( 7) 
0.2658( 7) 
0.2592( 8) 


0.2785( 7) 
0.2732( 7) 
0.2680( 8) 
0.2628( 8) 


0.4279(27) 
0.4240(28) 
0.4202(30) 
0.4165(34) 


0.3041(10) 
0.2993(10) 
0.2945(11) 
0.2899(12) 


0.3537(28) 
0.3482(30) 
0.3428(32) 
0.3374(35) 


0.3537(28) 
0.3487(30) 
0.3441(34) 
0.3406(40) 


0.3460(32) 
0.3407(33) 
0.3355(36) 
0.3306(40) 


0.3771(36) 
0.3728(39) 
0.3690(43) 
0.3661(48) 


0.1235 
0.1235 
0.1235 
0.1235 


0.1230 
0.1235 
0.1240 
0.1245 


0.2675( 7) 
0.2610( 8) 
0.2545( 8) 
0.2477( 8) 


0.2664( 7) 
0.2610( 8) 
0.2557( 8) 
0.2504( 9) 


0.4203(30) 
0.4168(33) 
0.4133(36) 
0.4099(40) 


0.2944(11) 
0.2896(12) 
0.2848(12) 
0.2801(14) 


0.3434(33) 
0.3378(35) 
0.3322(38) 
0.3268(42) 


0.3430(32) 
0.3378(35) 
0.3330(39) 
0.3291(46) 


0.3354(35) 
0.3299(38) 
0.3246(41) 
0.3195(45) 


0.3686(42) 
0.3643(45) 
0.3604(49) 
0.3575(56) 


0.1240 
0.1240 
0.1240 
0.1240 


0.1230 
0.1235 
0.1240 
0.1245 


0.2564( 8) 
0.2498( 8) 
0.2430( 9) 
0.2359( 9) 


0.2538( 8) 
0.2484( 8) 
0.2430( 9) 
0.2375(10) 


0.4129(36) 
0.4099(40) 
0.4070(44) 
0.4041(50) 


0.2849(13) 
0.2800(13) 
0.2751(14) 
0.2704(16) 


0.3332(40) 
0.3275(42) 
0.3218(46) 
0.3162(52) 


0.3324(38) 
0.3270(41) 
0.3218(46) 
0.3175(55) 


0.3247(41) 
0.3192(44) 
0.3136(48) 
0.3083(54) 


0.3608(50) 
0.3565(54) 
0.3527(60) 
0.3498(68) 


0.1245 
0.1245 
0.1245 
0.1245 


0.1230 
0.1235 
0.1240 
0.1245 


0.2451( 9) 
0.2383(10) 
0.2312(10) 
0.2237(11) 


0.2404( 9) 
0.2349( 9) 
0.2294(10) 
0.2237(11) 


0.4054(46) 
0.4035(51) 
0.4016(58) 
0.4000(68) 


0.2754(15) 
0.2705(16) 
0.2656(18) 
0.2609(19) 


0.3239(53) 
0.3180(57) 
0.3122(63) 
0.3065(71) 


0.3225(49) 
0.3168(53) 
0.3112(60) 
0.3065(71) 


0.3143(52) 
0.3085(56) 
0.3028(61) 
0.2972(70) 


0.3548(64) 
0.3507(69) 
0.3470(77) 
0.3444(89) 


fit 


range 


28-56 


28-56 


16-24 


28-58 


26-44 


26-44 


26-44 


26-44 



TABLE VIL The same resuhs as Table El for P = 6.10. 











Positive-parity baryons 






Negative-parity baryons 








"^oct(E) 


m-oct(A) 


TTlsing 




n^oct{E) 


n^oct{A) 




rridcc 


0.1230 


0.1230 


0.2382( 


4) 


0.2382( 


4) 


0.3466(30) 


0.2560( 7) 


0.2967(23) 


0.2967(23) 


0.2925(20) 


0.3103(26) 


0.1230 


0.1235 


0.2320( 


5) 


0.2329( 


5) 


0.3417(32) 


0.2511(08) 


0.2911(25) 


0.2919(25) 


0.2875(21) 


0.3059(28) 


0.1230 


0.1240 


0.2257( 


5) 


0.2276( 


5) 


0.3365(35) 


0.2463(08) 


0.2854(27) 


0.2874(27) 


0.2827(22) 


0.3017(32) 


0.1230 


0.1245 


0.2193( 


5) 


0.2223( 


5) 


0.3311(39) 


0.2417(09) 


0.2796(29) 


0.2835(33) 


0.2783(25) 


0.2981(38) 


0.1235 


0.1230 


0.2271( 


5) 


0.2262( 


5) 


0.3370(35) 


0.2462( 8) 


0.2867(26) 


0.2858(26) 


0.2826(22) 


0.3014(31) 


0.1235 


0.1235 


0.2208( 


5) 


0.2208( 


5) 


0.3325(38) 


0.2413(09) 


0.2809(28) 


0.2809(28) 


0.2776(24) 


0.2969(34) 


0.1235 


0.1240 


0.2143( 


5) 


0.2154( 


5) 


0.3277(42) 


0.2365(10) 


0.2750(31) 


0.2762(31) 


0.2728(26) 


0.2926(38) 


0.1235 


0.1245 


0.2076( 


6) 


0.2101( 


6) 


0.3226(48) 


0.2319(11) 


0.2689(34) 


0.2722(38) 


0.2683(29) 


0.2889(46) 


0.1240 


0.1230 


0.2160( 


5) 


0.2137( 


5) 


0.3268(42) 


0.2366(10) 


0.2768(32) 


0.2746(31) 


0.2729(26) 


0.2929(39) 


0.1240 


0.1235 


0.2095( 


6) 


0.2082( 


5) 


0.3232(47) 


0.2317(11) 


0.2709(34) 


0.2695(34) 


0.2680(28) 


0.2883(43) 


0.1240 


0.1240 


0.2028( 


6) 


0.2028( 


6) 


0.3193(53) 


0.2270(12) 


0.2646(37) 


0.2646(37) 


0.2631(31) 


0.2840(49) 


0.1240 


0.1245 


0.1958( 


6) 


0.1974( 


6) 


0.3150(62) 


0.2225(13) 


0.2581(42) 


0.2604(45) 


0.2586(35) 


0.2803(59) 


0.1245 


0.1230 


0.2049( 


6) 


0.2006( 


6) 


0.3151(55) 


0.2276(12) 


0.2680(44) 


0.2628(42) 


0.2638(33) 


0.2855(56) 


0.1245 


0.1235 


0.1982( 


7) 


0.1950( 


6) 


0.3133(63) 


0.2228(14) 


0.2617(48) 


0.2575(45) 


0.2588(36) 


0.2809(63) 


0.1245 


0.1240 


0.1912( 


7) 


0.1894( 


7) 


0.3114(74) 


0.2182(15) 


0.2551(53) 


0.2523(50) 


0.2540(40) 


0.2765(72) 


0.1245 


0.1245 


0.1839( 


8) 


0.1839( 


8) 


0.3094(91) 


0.2139(18) 


0.2479(61) 


0.2479(61) 


0.2493(46) 


0..2727(88) 


fit 


range 


40-64 


40-64 


24-36 


44-64 


36-52 


36-52 


32-52 


36-52 



the same size (~ 2 fm), and we cannot exam- 
ine the finite volume effects on these lattices. In 
Ref. 0, the finite volume effect on the negative- 
parity baryon masses was evaluated as 5% by com- 
paring the masses on lattices with volume sizes 
1.5fm and 2.2fm (1.6fm and 2.1fm) at /3 = 6.0 (6.2) 
on quenched isotropic lattices. This amount of fi- 



nite size effect may also exist in our results, while 
our lattice volumes are close to the larger ones in 
Ref. 0]. 

• Lattice discretization error 

Table lVIIll shows the baryon masses on each lattice. 
We find it hard to take the continuum limit even 
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TABLE VIII: The Hadron spectrum expressed in the unit of 
GeV at the physical quark masses. The K* meson mass is 
used for the determination of the scale unit a-r. 



0.4 0.6 
<mp,^(m> [GeV'l 

FIG. 2: The spectrum of the vector meson plotted against 
the pseudoscalar meson mass squared mps in the physical 
unit. The open symbols denote the direct lattice data, and 
the filled symbols denote the results for the physical quark 
masses obtained from the lattice data with the linear chiral 
extrapolation. 



for the ground-state baryons and mesons, since only 
the three /?'s are taken here and their behavior is 
not so smooth that one can apply a simple extrap- 
olation. In addition, such fluctuating behavior of 
data may be large due to the lack of statistics, gen- 
uine discretization errors would not be negligible. 
In the range of lattice cutoff 1-2 GeV, the fluctu- 
ation of masses is at most about 5%, except for 
the case of A^^^^ for which 7% deviation is found. 
This gives us a hint on the potential size of the 
discretization errors. We also note that Ref. 
examined how the 0{aa) and 0{a?) discretization 
effects decrease as /3 increases in the meson sectors, 
and found that those in the calibration of 7f are 
sufficiently reduced already ai (3 — 6.10. For these 
reason we discuss physical consequences of our re- 
sult mainly based on the data of /3 = 6.10 lattice 
in the next section. 

• Chiral extrapolation 

From the study of the chiral perturbation theory, 
there appears a non-analyticity in the chiral extrap- 





/3 = 5.75 


P = 5.95 


P = 6.10 


p 


0.7973(52) 


0.7965(53) 


0.8005(65) 


K* 


0.8939 


0.8939 


0.8939 


<P 


0.9905(55) 


0.9913(53) 


0.9873(66) 


N 


1.1118(84) 


1.0781(81) 


1.1055(72) 


A 


1.2095(75) 


1.1825(75) 


1.2002(67) 


E 


1.2256(76) 


1.1982(76) 


1.2173(67) 




1.3393(71) 


1.3183(75) 


1.3291(67) 


A 


1.366(23) 


1.342(16) 


1.3685(17) 


E* 


1.459(20) 


1.440(14) 


1.4586(15) 




1.552(17) 


1.538(12) 


1.5486(13) 


n 


1.645(15) 


1.635(11) 


1.6387(12) 


iV{-) 


1.686(79) 


1.599(59) 


1.618(57) 




1.784(67) 


1.705(49) 


1.717(49) 


=;(-) 


1.882(56) 


1.810(40) 


1.816(42) 




1.788(66) 


1.703(48) 


1.700(49) 


A'r) 


1.732(28) 


1.646(49) 


1.725(39) 




1.806(84) 


1.877(73) 


1.833(80) 




1.896(72) 


1.955(61) 


1.913(69) 




1.986(60) 


2.032(50) 


1.994(57) 




2.077(49) 


2.109(39) 


2.074(46) 


a(+' 

smg 


2.288(58) 


2.292(46) 


2.150(70) 



olation at the quenched level near the chiral limit 
as < 400 MeV 27]. We have, however, taken 
the naive linear extrapolation for both the positive- 
and negative-parity baryons, because the results for 
the baryon masses seem to be well described with 
the linear relation in the quark-mass region corre- 
sponding to ttItt > 600 MeV in the present calcu- 
lation. In addition, the behavior of the negative- 
parity baryon masses near the chiral limit is less 
known. Therefore, it is difflcult to estimate the 
non-analyticity in the chiral extrapolation from the 
present results. The quantitative estimate of its ef- 
fect on the negative-parity baryon masses is a fu- 
ture problem with high precision data with small 
quark masses. 

Quenching effects 

There is about 10% uncertainty for the ground- 
state hadron spectra coming from the quenching 
effect. Also for the negative-parity baryons, there 
should appear such a quenching effect. In addition, 
there may appear nontrivial excitation effect of the 
rj' meson in quenched QCD, where rj' degenerates 
with the other Nambu-Goldstone bosons in the chi- 
ral limit due to the ignorance of the fermionic de- 
terminant. Such an effect from r/ is reported to 
appear near the chiral limit as < SOOMeV [28|. 
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o octet{A)(- 
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0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 
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CO 

E 



1.6 
1.4 
1.2 
1.0 



1 1 1 1 1 1 1 1 

M-10 




□ octet(+) 
A decuplet(+) 
o octet{-) 
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iA*-' 




i'-^A E '-' 
1N<-' ^ 




E* j^-^ 




1 . 1 . 1 . 1 





0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0 
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□ octet(A)(+) 
A singlet(+) 
o octet(A)(- 
V singlet(-) 



0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0 
<mps'(mi)> [GeV'] 



FIG. 3: The spectra of positive- and negative-parity baryons plotted against the pseudoscalar meson mass squared m^s- For 

each /3, the octet and the dccuplct baryons arc shown in the left panel and octct(A) and singlet baryons in the right panel. The 
horizontal axis denotes the averaged pseudoscalar mass square. The open symbols denote the direct lattice data, and the filled 
symbols denote the results for the physical quark masses obtained from the lattice data with the linear chiral extrapolation. 
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IV. DISCUSSION 

Our numerical results for the hadron spectra are sum- 
marized in Fig. El and Table IVinl The masses of the 
negative-parity baryons are found to be heavier than 
those of the corresponding positive-parity sectors, as ex- 
pected. The flavor-singlet baryon is, however, an excep- 
tion: the positive-parity baryon is much heavier than 
the negative-parity one. This tendency seems consis- 
tent with the 3Q state in the quark model, in which the 
flavor-singlet positive-parity baryon belongs to the 70- 
dimensional representation of the SU{6) symmetry with 
the principal quantum number N = 2. This multiplet is 
in general heavier than that belonging to the negative- 
parity baryons. QCD sum rule analysis jl^ and the other 
recent lattice calculation also predict the mass of the 
flavor-singlet negative-parity baryon lighter than that of 
the positive one. 

In order to compare our lattice results with experimen- 
tal values, various baryon masses a,t (3 — 6.10 together 
with the experimental values are shown in Fig.^ For the 
positive-parity baryons, the nucleon and the delta masses 
are somewhat higher than the experimental ones. Note 
again that quenched QCD exhibits the non-analytic be- 
havior in the chiral extrapolation on the nucleon near the 
chiral limit of 111,, < 400 MeV In comparison with 
the naive linear extrapolation, this effect lowers the nu- 
cleon mass in the chiral limit, although we have not taken 
into account the non-analytic behavior because of the ab- 
sence of its signal in our relatively heavy quark-mass re- 
gion of ttItt > 600MeV. The other positive-parity baryons 
with strangeness reproduce the experimentally observed 
masses within 10% deviations. The better reproduction 
of strange baryon masses may be natural because of the 
following reason. For the strange baryon, the ambiguity 
from the chiral extrapolation would be less than that of 
the nucleon and the delta, because the strange quark is 
relatively massive and the non-analytic behavior arises 
only from up and down quark mass region. 

As for the negative-parity baryons, most of the present 
lattice results comparatively well reproduce the experi- 
mental spectra as shown in Fig. ^ in spite of a rela- 
tively large statistical error. However, the flavor-singlet 
negative-parity baryon is exceptional, and its calculated 
mass of about 1.7GeV is much heavier than the exper- 
imentally observed A(1405) with the difference of more 
than 300 MeV. The difference between the lattice result 
of 1.7GeV and the experimental value of the A(1405) is, 
however, the largest in all the hadrons in consideration. 
Even taking the quenching effect into account, this dis- 
crepancy seems too large. (Note again that the flavor- 
singlet baryon has one strange quark, and the ambiguity 
from the chiral extrapolation is expected to be less than 
that of the nucleon and the delta.) 

If the A(1405) could be well described as a three va- 
lence quark system, its mass would be reproduced with 
the simple three-quark operator in quenched QCD. How- 
ever, the light mass of the A(1405) is not reproduced 



in the present simulation. Therefore, the present lat- 
tice QCD result physically indicates that the experimen- 
tally observed A(1405) cannot be described with a simple 
three valence quark picture, that is, the overlap of the 
A(1405) with the 3Q component is rather small. This 
seems to support other possible pictures for the A(1405) 
such as the penta-quark state or the NK molecule. In 
this sense, lattice QCD simulations with the penta-quark 
operator would be meaningful to elucidate the nature of 
the A(1405). 

Here, we add several comments and cautions in 
quenched QCD. First, there is possible mixing between 
3Q and 5Q states through the "Z-graph" even in the 
quenched approximation 30]. Therefore, to be strict, the 
separability into 3Q and 5Q states does not hold even in 
the quenched approximation, although the mixing be- 
tween the two states is rather suppressed. Nevertheless, 
the conclusion of the non-3Q picture for the A(1405) is 
still plausible, because, if the A(1405) is described as the 
3Q state, its mass is to be reproduced with the 3Q oper- 
ator in quenched QCD. Second, for the definite conclu- 
sion, we have to pay attention to the non-analytic be- 
havior in the flavor-singlet negative-parity baryon near 
the chiral limit (8|, in spite of the naive expectation of 
its smaller effect for strange baryons. Third, according 
to the neglect of the fermionic determinant, 77' becomes 
unphysically "light" as the Nambu-Goldstone particle in 
quenched QCD, and the non-unitary behavior appears 
in the baryon correlator due to the "light" ry' excita- 
tion near the chiral limit below ^ 250 MeV [28l |. 
where there appears an unphysical "decay" process of 
the negative-parity baryon into an t^'-N state. Although 
these non-analytic and the non-unitary behaviors are not 
observed in the present simulation with relatively heavy 
quark masses as rriTr > 600MeV, these effects should be 
taken into account for the simulation near the chiral limit. 

We now focus on other negative-parity baryons. The 
mass ratio between the positive- and the negative-parity 
baryons is shown in Table IIXI For both the octet and 
the decuplet baryons, the relative mass difference be- 
tween the parity partners becomes smaller, as the aver- 
aged quark mass increases by the inclusion of the strange 
quark. This tendency is experimentally observed for the 
octet baryons. (The empirical identification of negative- 
parity decuplet baryons is not established.) This behav- 
ior is also reported in another lattice QCD analysis by 
the domain wall fermion [3. From Fig, m and Table Hxl 
we find that the lattice results of the flavor octet and de- 
cuplet baryons are all close to the observed lowest-lying 
negative-parity baryons, the A^(1535), A(1670), 11(1620) 
and A (1700), in spite of the relatively large statistical er- 
ror. The I](1620), which is experimentally confirmed as 
the negative-parity strange baryon with — l/2~ [T6| . 
is consistent with the parity partner of the E baryon. 
The parity partner of the S baryon is expected to be the 
S(1690) from our calculation, although the spin-parity 
of the S(1690) is not yet confirmed experimentally. Re- 
cently it has been proposed based on the chiral unitary 
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approach 'si'l that the S(1620) has the negative-parity, 
although its experimental status is still one-star (evi- 
dence of existence is poor). It is, however, difficult to 
distinguish between the S(1620) and the S(1690) from 
our result due to the statistical errors. For the decu- 
plet baryons, we can regard the parity partner of the 
A(1232) as the A(1700), although the experimental data 
is poor. The positive-parity flavor-singlet baryon is found 
to be much heavier than the negative-parity decuplet, 
and hence its investigation seems much difficult both the- 
oretically and experimentally. 

Finally, we comment on recent lattice studies on the 
negative-parity baryons. Sasaki et al. investigated 
the negative-parity non-strange baryon N^^\ the par- 
ity partner of the nucleon N^~^\ with the domain wall 
fermion 0. Their lattice is 163x32 at /3 = 6.0 {a'^ = 1.9 
GeV) and the result is N'-^^ /N^+^ ~ 1.45, which is con- 
sistent with ours. Gockler et al. studied the negative- 
parity non-strange baryon N^^^ using the 0(a) improved 
Wilson quark on the isotropic lattices with the size of 
16^ X 32 and 32^ x 64 @. They obtained the similar re- 
sult, N^-yN^+'f = 1.50(3). Melnitchouk el al. also stud- 
ied the negative-parity baryons using the 0(a) improved 
Wilson quark on the isotropic lattice, 16'^ x 32(a = 0.125 
fm) 0. Since they did not carry out the chiral extrap- 
olation, we do not compare the results quantitatively, 
but the qualitative behavior is similar to us. They also 
investigated the flavor-singlet baryons. Instead of the 
flavor-singlet interpolating field, they used the "common" 
interpolating field which is the common part of the in- 
terpolating fields for the octet A hyperon and the singlet 
baryon. The result is much heavier than the experimen- 
tal value of the A(1405) even for such an field. Lee et al. 
investigated the excited state baryons with the overlap 
fermion with the lattice 16'^ x 28 They employed the 
constrained curve fitting method for the mass fitting and 
obtained the baryon masses lower than those from the 
conventional fitting method. Thus, their results seem 
to be lower than ours and the other's, while they did 
not carry out the chiral extrapolation. Dynamical quark 
simulation of excited-state baryons is also in progress . 
As for the negative-parity nucleon, their present result 
is consistent with the quenched result within statistical 
errors. 



V. SUMMARY AND CONCLUDING REMARKS 

We have studied the mass spectra of the negative- 
parity baryons and the flavor-singlet baryons in quenched 
anisotropic lattice QCD. We have used three lattices of 
almost the same physical spatial volume of about (2fm)3 
with the spatial cutoffs = 1-2 GeV and the renor- 
malized anisotropy ~ ag- / ar — -i. We have adopted the 
standard Wilson plaquette gauge action and the 0{a) 
improved Wilson quark action at the tadpole-improved 
tree-level j2j|. The positive- and negative-parity baryon 
masses are extracted with the parity projection from the 



same baryon correlators based on the three valence quark 
picture. 

For the flavor octet and decuplet negative-parity 
baryons, the calculated masses are close to the exper- 
imental values of corresponding lowest-lying negative- 
parity baryons. For several negative-parity baryons, our 
lattice data have suggested some predictions. For in- 
stance, A(1700) can be regarded as the parity partner of 
A(1232), and the S(1690) would be the parity partner of 
the 5 baryon, although the spin-parity of the S(1690) is 
not yet confirmed experimentally. 

As for the flavor-singlet negative-parity baryon, such a 
three-quark state has been found to lie around 1.7GeV, 
and has been much heavier than the A(1405). Even con- 
sidering the systematic errors, this difference of about 
300 MeV seems too large. If the A(1405) is described 
as a three valence-quark state, its mass would be repro- 
duced in the present simulation. In fact, the present lat- 
tice result which cannot reproduce the A(1405) physically 
implies that the A(1405) is not described as the simple 
three quark picture, i.e., the overlap of the A(1405) with 
the three-quark state is rather small. This seems to sup- 
port an interesting picture of the penta-quark state or the 
NK molecule for the A(1405). For more definite under- 
standing of the A(1405), it would be desired to perform 
lattice QCD simulations in terms of the NK molecule or 
the penta-quark state. Such a study is interesting even at 
the quenched level, where dynamical quark loop effect is 
absent and then the quark-level constitution of hadrons is 
clearer. As for the positive-parity flavor-singlet baryon, 
its calculated result is found to be much heavier than 
the negative-parity one, as is consistent with t he q uark 
models [iflllll and the QCD sum rule analysis |29|. 

Very recently, LEPS Collaboration has experimentally 
observed the 0+ (or Z~^) baryon with 5 = -|-1 [23|j which 
requires at least "five valence quarks" as uudds and is 
physically identified as a "penta-quark system" . The 
comparison between the baryon and the A(1405) may 
be useful to investigate the features of the penta-quark 
system. It is also interesting to investigate this type of a 
penta-quark system using lattice QCD simulation. 

From aspect of the chiral symmetry, the parity partner 
should be degenerate if the symmetry is restored at finite 
temperature and/or density. It is interesting to see how 
the mass difference between the parity partners changes 
at finite temperature on the lattice. Several works on it 
are already reported for the screening mass of the nu- 
cleon [s^ and they favor the parity degeneracy at the 
chiral phase transition. Recently based on the chiral ef- 
fective theory such as the linear sigma model [34] and the 
chiral perturbation theory , two different assignments 
for the negative-parity baryons have been proposed: un- 
der the chiral transformation, the negative-parity baryon 
transforms in the same way as the positive-parity one 
in one scheme and in the opposite way in the other. 
These two assignments behave differently toward the chi- 
ral restoration. Therefore, it is interesting to study them 
from the quark degrees of freedom such as in lattice QCD 
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FIG. 4: Various baryon masses obtained from the /3 = 6.10 lattice. For the negative-parity baryons, the experimental values 
of iV(1535), A(1670), E(1620), H(1690), A(1700) and A(1405) are added. 



TABLE IX: The ratio of negative and positive-parity baryon masses. In the last line, the ratio of flavor-singlet negative-parity 
and octet positive-parity baryon masses is also listed. Physical values of the negative-parity baryons are taken to be iV(1535), 
E(1620), H(1690), Aoct(1670), Asi„g(1405) and A(1700). (*)Note that the spin-parity of the H(1690) is not yet confirmed. 





f3 = 5.75 


13 = 5.95 


13 = 6.10 


physical value 


Ar(-)/iv(+) 


1.516(70) 


1.484(54) 


1.463(51) 


1.635 


S(-)/E(+) 


1.456(54) 


1.423(40) 


1.410(40) 


1.358 


=(-)/=+ 


1.405(41) 


1.373(30) 


1.366(31) 


1.282* 


a(-) /a(+) 

-"^oct / -"-Oct 


1.479(54) 


1.440(40) 


1.417(40) 


1.496 


a'.-) /a(+) 


0.757(22) 


0.718(26) 


0.802(32) 




A(-)/A(+) 


1.322(64) 


1.399(55) 


1.339(61) 


1.380 


s*H/s*(+) 


1.299(52) 


1.358(43) 


1.312(49) 




H*H/H*(+) 


1.280(41) 


1.321(33) 


1.288(38) 






1.263(31) 


1.290(24) 


1.266(29) 




smg oct 


1.432(23) 


1.392(42) 


1.437(34) 


1.259 



at finite temperature. 
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